The EPR experiment in the energy-based stochastic 
reduction framework 



J Silman 1 , S Machnes 1 , S Shnider 2 , L P Horwitz 1 3 and A 
Belenkiy 2 

1 School of Physics and Astronomy, Tel- Aviv University, Tel- Aviv, Israel 
2 Department of Mathematics, Bar-Ilan University, Ramat-Gan, Israel 
3 Department of Physics, Ariel University Center of Samaria, Ariel, Israel 

Abstract. We consider the EPR experiment in the energy-based 
stochastic reduction framework. A gedanken set up is constructed to 
model the interaction of the particles with the measurement devices. 
The evolution of particles' density matrix is analytically derived. We 
compute the dependence of the disentanglement rate on the parameters 
of the model, and study the dependence of the outcome probabilities 
on the noise trajectories. Finally, we argue that these trajectories can 
be regarded as non-local hidden variables. 



1. Introduction 

A pure quantum state of a system is a vector in a Hilbert space, which may be 
represented as a linear combination of a basis of eigenstates of an observable (self- 
adjoint operator) or of several commuting observables. Let us suppose that the 
eigenvalues corresponding to the eigenstates of the Hamiltonian operator of a system 
are the physical quantities measured in an experiment. If the action of the experiment 
is modeled by a dynamical interaction induced by a term in the Hamiltonian of the 
system, and its effect is computed by means of the standard evolution according to the 
Schrodinger equation, the final state would retain the structure of the original linear 
superposition. One observes, however, that the experiment provides a final state that 
is one of the basis eigenstates and the superposition has been destroyed. The resulting 
process is called reduction or collapse of the wave function. 

The history of attempts to find a systematic framework for the description of this 
process goes back very far in the development of quantum theory (e.g., the problem of 
Schrodinger's cat p]). In recent years significant progress has been made. Rather than 
invoking some random interaction with the environment and attributing the observed 
decoherence, i.e. collapse of a linear superposition, to the onset of some uncontrollable 
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phase relation, more rigorous methods have been developed. These methods add to 
the Schrodinger equation stochastic terms corresponding to Brownian fluctuations of 
the wave function [21 [3], [5l El [Tj [HI [9| , generally understood as arising from the 
presence of the measurement device. 

In this paper, we apply some of these state reduction methods to the phenomena 
considered in Bohm's formulation [10] of the Einstein-Podolsky-Rosen paradox [TT] 
(henceforth the EPRB paradox), later analyzed by Bell for its profound implications 
|12| . and explored experimentally by Aspect et al. [13]. The system to be studied 
consists of a pair of spin-| particles in the singlet state 

liM = ^(IU>-UT», (i) 

the arrows denoting the spin components of the particles relative to some arbitrary 
axis. The determination of the spin state of one of the particles implies with certainty 
the spin state of the other, even when the particles are very far apart. The particles 
are therefore said to be entangled. 

The question is often raised as to how the state of the second particle can respond 
to the arbitrary choice of direction in the measurement of the first. This question is 
dealt with here by the construction of a gedanken set up describing the interaction 
of the particles with the measurement devices. On this basis, using the mathemat- 
ical models recently developed for describing the reduction, or collapse, of the wave 
function, we answer this question and give a mathematical description of the process 
underlying such a measurement. 

The paper is organized as follows. We begin in section 2 by reviewing the 
energy-based stochastic extension of the one-particle Schrodinger equation, and discuss 
its generalization to noninteracting multiparticle systems. In section 3 we present 
a gedanken set up for studying the EPR experiment and show that it leads to 
the expected quantum mechanical predictions. Next, in section 4 we analytically 
compute the stochastic expectation of the particles' density matrix, and quantify 
their disentanglement rate. In section 5 we simulate the evolution of the state of 
the particles for different random realizations of the noise, and argue that the noise 
trajectories can be regarded as nonlocal hidden- variables. We end by discussing future 
avenues of research. 
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2.1. Energy-based stochastic extension of the Schrddinger equation 

In the energy-based stochastic reduction framework the Schrodinger equation is 
extended as follows [3 [9] 

d\m) = -iH\m)dt-^{H-H(t)y\m)dt 

+ "-{H-H{t))m))dW(t). (2) 

Here H(t)={ip(t)\H\ip(t)) (we assume throughout that \i^(t)) is normalized, in 
consistence with eq. ([2|) as it is norm preserving, W(t) is a standard Wiener process, 
and ^ is a parameter characterizing the reduction time scale the rate. (Note the choice 
of "natural" units h = c = 1. Accordingly, all quantities throughout the paper are 
expressed in units of length [£].) 

From the Ito calculus rules it immediately follows that the above process has two 
basic properties 

(i) Conservation of energy 

H(t) = H(0)+<; [ dW(s)V(s), (3) 
Jo 

where V(t)=(ip(t)\(H — H(t)) 2 \ip(t)) is the variance of the energy process H(t). 

(ii) Stochastic reduction 

dV(t) = -<, 2 V 2 {t)dt + <;/3(t)dW(t) , (4) 
where (3(t)=(ip(t)\(H — H(t)) \tp(t)) is the third moment of the energy deviation. 

It follows from eq. |(4]) that the expectation E[V(t)) of the variance process obeys the 
relation [3[8] 

E[V(t)] < E[V(0)] - e [ dsE[V(s)} 2 . (5) 
Jo 

Since V(t) is positive, this implies that E[V(t — » oo)] — > and (up to measure zero 
fluctuations) V(t -> oo) -> 0. And since V{t) = {ip{t)\{H - H{t)) 2 \i){t)) , V(t) = 
implies (ip(t)\(H - H(t))\<ip(t)) = or H\ip(t)) = H(t)\tp(t)), so that \tp(t)) is an 
eigenstate of the Hamiltonian. Assuming no degeneracy, the system therefore reduces 
to one or another of the eigenstates of the Hamiltonian H, in accordance with the 
statistical predictions of standard quantum mechanics [8], Therefore, the expectation 
of the final configuration E[\ip(t — * oo)) (ip(t — > oo)|] corresponds to a mixed state, 
with each of the diagonal elements an eigenstate of H. 

Note that the framework we have described cannot differentiate between 
degenerate eigenstates. When this is the case, as in the standard theory [15] . the 
reduction process drives the system to degenerate subspaces with the original relative 
phase between the spanning eigenstates remaining unchanged. 
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2.2. Extension to noninteracting multiparticle systems 

Nothing in the previous subsection limits the discussion to single particle systems. 
The Hamiltonian in eq. ([2]) may just as well represent a multiparticle system. This, 
however, is not the only possible generalization to multiparticle systems, and indeed 
there are cases where it is not suitable. To see this, and in anticipation of the next 
section, let us consider a pair of noninteracting particles A and B. The Hamiltonian 
is now a direct sum 

H = H A ®H B . (6) 

We assume that the particles are very far apart, and that the environment does not 
carry pervasive long-range correlations. Under these conditions the evolution of the 
state's stochastic expectation E[p(t)] {pit) = \ip(t)) should be local, in the sense 

that no correlations, quantum or classical, are generated. 

Bearing this in mind, let us plug the Hamiltonian eq. ((6|) into eq. |(2]). Averaging 
over the noise we obtain a Lindblad type equation [16} [17] for the state's stochastic 
expectation 

±E[p(t)]=-i[H,E[p(t)}}-'j J2 [Hi,[Hj,E[p(t)]]]. (7) 

i,j=A,B 

This equation is causal (does not allow for superluminal signalling), as is easily 
established by tracing over any of the two subsystems. However, the mutual 
information^] may increase with time. The evolution is therefore nonlocal, as may 
well have been expected considering that both systems are driven by the same noise. 

However, a local evolution equation for the state's stochastic expectation can be 
achieved if we have each of the systems driven by an independent noise term. This 
means that eq. J2]) must be generalized as 

d\i/>(t)) = -iH\i/>(t))dt-^ J2 ^(Hi-Hi^fl^dt 

i=A, B 

+ \Y. -#<(*)) hK*)><Wi(*) (8) 

i=A, B 

where q governs the reduction rate of particle i to the eigenstates of Hi and 
dWi(t)dWj(t) — Sijdt. Indeed, the above process invariably drives the system to 
product states of the form \Ea) <8> \Eb), where ff, \Ej) = E { \Ej), with the same 
probabilities as predicted by the standard theory [8]. The corresponding Lindblad 
equation for the state's stochastic expectation is now given by 

jEW)\ = ~i[H, E[p(t)]} - [ ^' E[pm ■ (9) 

2—1, 2 

I The mutual information of two systems is defined as Iab = Sa + Sg — Sab j where Si and Sab 
are the von-Neumann entropies of system i and the composite system, respectively, and serves as a 
quantitative measure of the total amount of correlations, quantum and classical, between the systems. 
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Note that in standard quantum theory this evolution can only arise from the respective 
coupling of a pair separate noninteracting systems to noncorrelated environments. 

3. Stochastic reduction in the EPRB experiment 

In this section we construct a gedanken set up to show how the energy-based stochastic 
reduction framework can provide us with a consistent description of the EPRB 
experiment. 

We consider a pair of spin-half particles, A and B, with vanishing total spin 
and momentum, moving in opposite directions. Along the path of each particle a 
spin measurement device is placed. With no loss of generality we assume that the 
measurement device in the path of particle A measures its spin component along z, 
and that the measurement device in the path of particle B measures its spin component 
along n, where h is some unit vector pointing in an arbitrary direction. 

So long as the particles are far from the measurement devices the Hamiltonian 
governing their (free) evolution is given by 

2mA 2rriB 

However, once the particles approach some neighborhood of the detectors (which we 
assume happens simultaneously) we assume that Hq is corrected by the addition of a 
perturbation 

Hm = Ha0 a ® 1 + Maimer! (11) 

describing the local interaction of the particles with the measurement devices. Here fa 
denotes the strength of the coupling of particle i to the corresponding measurement 
device and a B = er ■ n. The eigenstates of the perturbed Hamiltonian H = H + H int 
are products of momentum and spin eigenstates, and are fully specified by the four 
eigenvalues — oo < pa, Pb < oo and a A , a B — ±1. 

The continuous spectrum of the momentum operators gives rise to an irremovable 
degeneracy in H. Nevertheless, for wave packets localized in momentum space and 
sufficiently large values of the fa this residual degeneracy is negligible. In all that 
follows we shall therefore make the approximation 

H ~ H int ■ (12) 

The four possible spin outcomes are then given by It/ 1 ), It,/), II/), and ||/), the 
slanted (vertical) up and down arrows denoting spin-up and spin-down eigenstates of 
a B ( a A)> respectively, with corresponding probabilities 



1 2 9 1 o 9 

p \/ = p iy = g cos 2 ' T/ = Ly = 2 Sm 2 ' 



(13) 



The EPR experiment in the energy-based stochastic reduction framework 



6 



where 6 is the angle between the z and h axes. The reduction process ultimately 
reproduces the results of the standard theory. 

It is important, lending further credibility to our set up, that the measurement 
of the spin of only one of the particles suffices to induce the reduction. This is clearly 
seen by recasting H int in its diagonal form 

Hint -EE ^ ly> fa'l > ( 14 ) 

where A-|y = —X^y = )ma + Pb and \\y = —\[y = pa — Pb- Indeed, this is 
approximately what will be observed at a time t ~ 1/q when q -C . 

We also note that when pa equals Hb ( — Ms) the degeneracy of the states |tvO 
and |J./*) (IT/ 1 ) an d II./)) is not removed. This, however, is merely an artifact of not 
taking into account the quantum nature of the measurement devices and the fields 
generating the coupling. Indeed, the local uncertainty of these fields, together with 
the absence of significant long-range correlations, serves to lift this degeneracy. 

Finally, we should remark that if the two particles are identical, then the 
Hamiltonian must be symmetric under the interchange of the particle indices. 
However, since the particles are very far apart when the measurement takes place, 
there is no overlap of the one particle wave functions, and the symmetrization or anti- 
symmetrization of the composite wave function is not required. Thus, the presence 
of two widely separated measurement devices can split the degeneracy into distinct 
states, which can, in fact imply that Hint is not symmetric under particle exchange. 

4. Evolution of the state's stochastic expectation 

To study the evolution of the state's stochastic expectation in our set up we simply 
substitute pa<J% an d Pbo"% for Ha and Hb in eq. ([9]). We thus have 

jE[p{t)] = -iHAlvl, E[p(t)]\ - ^[a%, [a% E[p(t)}}} 

- in B [cr%, E[p{t)]] - ^-{a%, [a%, E[p(t)}]] . (15) 

This equation is linear in E[p(t)\, and can therefore be transformed into a linear 
equation for a vector whose components are the sixteen elements of E[p(t)]. An 
analytical solution can then be obtained by bringing the 16 x 16 matrix, representing 
the action of the operators on the right-hand side of the equation, to its Jordan normal 
form. When working in the basis \y , I/ 1 , \y the solution reads as follows 





( sin2 1 
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Figure 1. (Color online) The negativity as a function of the 
coupling strengths /ja and /is. (?a = fs = 1 M 1//2 , t = 1 [(], 
and 6 = 3tt/5) On the (dark blue) flat surface the negativity 
equals zero. The disentanglement time is finite unless fiA or 
fiB vanishes. 

(16) 

with Qi=exp (pi (4i — PiS?) t/2) (i = A, _B), and 9 the angle between z and n. We 
note that the Qi can be decomposed into a product of exp (2i/ijt) and exp (— /if <^t/2) 
representing the contibutions of the unitary and stochastic processes, respectively. We 
also note that in the limit that t — - > oo the off-diagonal terms vanish and the density 
matrix reproduces the expected measurement outcomes and probabilities. 

We now wish to examine the transition of the stochastic expectation of the density 
matrix from the initial maximally entangled singlet state to the mixture of the final 
outcomes, paying particular attention to the rate of the disentanglement. While for 
pure states, the entanglement is quantified by the von-Neumann entropy, there is no 
single measure of entanglement for mixed states. One of the standard measures is the 
negativity |18j . 5(||p Ti || — 1), where p Ti , the partial transpose with respect to i of p, is 
obtained from p by transposing the indices of system |§[ which is just minus the sum 
of the negative eigenvalues of the p Ti [IHl [20] . 

Numerical results are presented in figures (1-3). Figure (1) illustrates the depen- 
dence of the negativity on the coupling strengths, pa and ps- We see that unless pa 
or pb vanishes, the disentanglement time is finite. This phenomenon, termed "entan- 
glement sudden death" [21) . is not unique to our setting and is typical of open systems 
dynamics [22J. For short times the rate of disentanglement is roughly dependent on 
\J pa 2 + Pb 2 , while for long times becomes linear. Figures (2) and (3) show the neg- 
ativity as function of time for different values of <,a = and 8, respectively. 

§ When considering the entanglement between two systems it is irrelevant which of the indices is 
transposed. 
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Figure 2. (Color online) 
The negativity as function 
of time for different values 
of the SLA = SB. {6 = 
3n/5, ha = MB = 
The legend gives the value of 
Si 2 /8 for each of the curves. 



Figure 3. (Color online) 
The negativity as function of 
time and 6. (sla = Sb = 
1M 1 / 2 , ha = Mb = 1) The 
legend gives the value of 9 
for each of the curves. 



5. Noise trajectories 

Eq. {3]) gives rise to an entropy conserving evolution. In particular, this means 
that a pure state remains pure. The trajectories realized by the Wi(t) during the 
evolution fully specify the state's history. However, as there is no means of determining 
these, all accessible information is contained in the state's stochastic expectation. 
The trajectories can therefore be regarded as hidden-variableljjl, i.e. indeterminable 
variables carrying information regarding the state of the system unavailable in the 
standard quantum mechanical description. 

Prom the relation dWi(t)dWj(t) = Sijdt it does not follow that the trajectories are 
local hidden- variables. Indeed, any hidden-variable theory adhering to the statistical 
predictions of standard quantum mechanics must violate some Bell inequality, and 
as such our set up is manifestly nonlocal. Explicitly, this just means that the both 
Wa (t) and Wb (t) determine the final state of each of the particles (rather than Wi (t) 
determining the final state of particle i on its own) . This point is illustrated in figures 
(4) and (5), which present the results of numerical iterative solutions to eq. (JSj) 
with randomly generated noise terms. Figure (4) explicitly shows how, for the same 
realization of Wa (t) , different realizations of Wb (t) lead to different outcomes in the 
spin measurement of particle A (and particle B). It is also interesting that in this 
case the probabilities for the measurement outcomes no longer agree with those of 

|| For a somewhat different approach see [14] . 
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Run 15. Collapse to "00" 



Run 16. Collapse to "01" 





Run 2. Collapse to "10" 



Run 19. Collapse to "11" 





Figure 4. (Color online) Four realizations of the collapse process. The 
probabilities for the four different outcomes (1= T, 0= |) are plotted as a function 
of time. n A = hb = ?a = ?s = 0.14 [e] 1 / 2 and 9 = 5tt/11. The same 

realization of Wjiit) was used in all of the runs, whereas Ws(i) varies from one 
run to the next. We see that the final outcome of particle A depends on not only 
W A (t) but on W B (t) as well. 



quantum mechanics, as is evident from figure (5). 

To see how this comes about we must go back to eq. ([8]). Even though it 
describes a pair of uncoupled systems, it gives rise to a potentially disentangling 
(and entangling) evolution, because each of the Hi{t) depends on the full state of the 
system (Hi(t) — trj^i[p(t)Hi\), and therefore on Wj^i(s < t). Note, however, that 
this nonlocality does not allow for superluminal signalling since the trajectories are 
"hidden". 



6. Some concluding remarks 

We have discussed Bohm's formulation of the EPR experiment in the the energy-based 
stochastic reduction framework. In particular, we have seen how the presence of the 
measurement devices induces the reduction of the singlet state to the expected outcome 
product states with correct probabilities as predicted by the standard theory and have 
given the explicit time evolution of the process of disentanglement. As an extension 
of this idea, one may consider a problem with a natural degeneracy of some initial 
state where the presence of effective detectors of some type induces a perturbation 
in which stochastic reduction takes place, as in the asymptotic cluster decomposition 
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Figure 5. (Color online) Dependence of the outcome probabilities on the noise 
trajectories. The first three figures display the results of 1000 runs simulations, 
in which Wa^) is held fixed from run to the next, but realized differently in each 
simulation, while Wsit) varies from one run to the next. We see that different 
realizations of Wa{1) lead to different outcome probabilities that do not agree 
with those of quantum mechanics. In the fourth figure, displaying the results of 
10000 runs simulation, both Wji(t) and Ws{t) vary, and the quantum mechanical 
predictions are obtained. 



of products of quantum fields reducing an n-body system to m fc-body systems, or 
the formation of local correlations in n-body systems such as liquids, or spontaneous 
symmetry breaking. In all these cases, due to the existence of continuous spectra, 
there will be some residual dispersion in the final state, although possibly very small. 
We are currently studying possible applications of the methods discussed here to such 
configurations. 
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